Robust Nodal Structure of Landau Level Wave Functions Revealed by Fourier 
Transform Scanning Tunneling Spectroscopy 
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Scanning tunneling spectroscopy is used to study the real-space local density of states (LDOS) of 
a two-dimensional electron system in magnetic field, in particular within higher Landau levels (LL). 
By Fourier transforming the LDOS, we find a set of n radial minima at fixed momenta for the nth 
LL. The momenta of the minima depend only on the inverse magnetic length. By comparison with 
analytical theory and numerical simulations, we attribute the minima to the nodes of the quantum 
cyclotron orbits, which decouple in Fourier representation from the random guiding center motion 
due to the disorder. Adequate Fourier filtering reveals the nodal structure in real space in some 
areas of the sample with relatively smooth potential disorder. 

PACS numbers: 73.43.-f, 73.22.-f, 73.43.Cd, 73.20.At 



Directly mapping the wave functions of electrons gives 
the most pertinent access to the quantum mechanical 
properties of matter (H-Q- Under a perpendicular mag- 
netic field B in two-dimensional electron systems (2DES), 
self-interference of the circular electronic orbits leads to a 
standing wave pattern of probability density, as first cal- 
culated by Landau Q . The kinetic energy becomes quan- 
tized into discrete Landau levels (LL) E n = Huj c (n + |), 
with uj c the cyclotron frequency, h Planck's constant, 
and n = 0, 1, 2 . . . characterizing the number of nodes 
in the LL wave functions. Experimental observation of 
this nodal structure has until now remained elusive. The 
2DES are usually deeply buried in semiconducting het- 
erostructures, which prevents the use of high resolution 
scanning tunneling spectroscopy (STS). Recently, the ad- 
vent of surface 2DES in doped semiconductors 0-0]? 
graphen e lid, lllj and on the surface of topological in- 
sulators [14|, [l5| has in principle opened the way to such 
direct high resolution measurements. This should allow 
probing the internal structure of the LL wave functions. 

However, the Landau energy levels are highly degener- 
ate, so that the associated wave functions will be strongly 
disturbed by any perturbation such as disorder. One has, 
thus, to deal with the inherent complexity of disorder in 
spectroscopic measurements. More fundamentally, dis- 
order is crucial for the understanding of universal quan- 
tized Hall conductance [l6| in 2DES. At high magnetic 
field, disorder essentially lifts the LL degeneracy keeping 
n as a good quantum number due to the large cyclotron 
gap. The electronic motion is then largely decomposed 



into independent fast cyclotron orbit and slow drifting 
motion of the guiding center along equipotential lines of 
the smooth disorder landscape. Most of the equipoten- 
tial lines in the bulk being closed, this picture provides 
a simple localization mechanism that ensures a reservoir 
of localized electronics states between LLs, ultimately re- 
sponsible for the formation of wide quantized plateaus of 
Hall conductance [6]. This semiclassical picture was re- 
cently confirmed by STS real space imaging of the elec- 
tronic probability density [H, . The latter was found 
in the lowest LL, LLO, to follow equipotential lines with 
a transverse spread on the scale of the magnetic length 
Ib = Within the higher LLs, LLn for n > 0, 

the drift motion is expected to be accompanied by larger 
cyclotron orbits characterized by the quantum Larmor 
radii R n = IbV^™ + 1- This larger spread of the mean- 
dering LL wave functions in the transverse direction to 
the guiding center motion should contain signatures of 
the LL nodal structure for n > 1. Note that these nodes 
can be related to phase singularities of LL wave func- 
tions with n a winding number [12|, [l3| . This topologi- 
cal origin confers some robustness to the nodal structure 
which should be viewed as a key property of quantum 
Hall states. The purpose of this Letter is to show how to 
reveal these interference effects. 

The quantum states in disordered LLs are delicate to 
analyze in detail, because drift and cyclotron motion are, 
strictly speaking, not disentangled. We will show that a 
useful spectroscopic analysis can still be made by two- 
dimensional (2D) Fourier transform of the STS data, 





FIG. 1. (color online) (a)-(g) Real-space LDOS taken at B = 6 T for the sample voltages V s marked as circles in (h); 
measurements are done by lock-in technique with the modulation voltage V m od — 1 mV after stabilizing the tip at a current 
^stab =0.1 nA and a sample voltage V^tab = 150 mV; all dl/dV images are recorded in the same area and are displayed using 
the same color scale; values of the chosen sample voltage are: V s = —117 mV (a), -103 mV (b), -79 mV (c), -73 mV (d), -57 
mV (e), -29.4 mV (f), -23.1 mV (g). White arrows at the same positions in (a), (c), (g) mark localized states, which exhibit 
additional nodal structure in (c), (g). Crosses in (b), (e), and (f) mark states localized at potential hills; dashed rectangle in 
(b) marks an area around a drift trajectory, (h) Spatially averaged dl/dV curve using the 1600 curves recorded in the area of 
(350 nm) 2 ; dashed lines: Gaussian fits to the two spin levels of LLO. 



which allows to deconvolute the contribution from the 
discrete LL orbits and from the disorder induced drift 
motion. The nodal structure of LLn is revealed by a 
succession of n ring-like patterns in momentum-space, 
given by a set of fixed and disorder-independent minima 
within the momentum scale R n /l 2 B . Analytical and nu- 
merical calculations are performed to vindicate our find- 
ings. Moreover, a methodology to show the nodal struc- 
ture in real-space is demonstrated. 

The 2DES was prepared by 1% monolayer Cs adsorp- 
tion on an n-type InSb(llO) surface [l7| with donor (ac- 
ceptor) density N D = 9-10 21 m~ 3 (N A = 5-10 21 m" 3 ) Q. 
The STS measurements [s[ were performed in ultra-high 
vacuum at B = 6 T and temperature T = 0.3 K [18]. Fig- 
ures QJa)-(g) show images of the differential conductivity 
dl/dV recorded within the same spatial area for seven 
different sample voltages V s spanning the four spin-split 
LLs from LLO to LL3 [see the spatially averaged d//dV 
curve shown in Fig. Hpi)]. These spatial maps repre- 
sent the local density of states (LDOS) consisting of all 
wave functions at energy E ~ eV s within the experi- 
mental energy resolution of 2.5 meV [!, [l9[, see discus- 
sion in [2(j. In the lower tail of LLO [Fig. QJa)], several 
spatially isolated closed loops are visible, which corre- 
spond to individual localized states encircling potential 
minima as shown in Ref. [8]. The same states are visi- 
ble within the lower tail of the higher spin level of LLO 
[Fig. OJb)], however, superimposed to the states marked 
by white crosses which belong to the upper tail of the 
lower spin level and encircle potential maxima. Similar 



patterns are found at the same positions in the lower tail 
of higher LLs [marked by arrows for LL1 in Fig. [TJc) 
and LL3 in Fig. Dig)], but they are spatially wider and 
are, thus, more difficult to discriminate. The widening 
of the drift states marks the progressive increase of R n 
with n. Careful inspection of the localized states in the 
LL's lower tail [Figs. QJa) , (c) , (g) , see arrows] and in the 
upper tail [Figs.QJb),(e),(f)] reveals that the closed loops 
exhibit more complex oscillatory patterns perpendicular 
to the loop in higher LLs. 

In order to relate the oscillatory features to the nodal 
structure of the LL wave functions, we firstly overcome 
the difficulty that the cyclotron motion is randomly cor- 
related to the guiding center motion. This partly masks 
the nodal structure within our real-space data. Note that 
LL mixing would also blur this structure. However, in the 
present experimental conditions (large cyclotron gap and 
smooth disorder, see Ref. LL mixing plays a negli- 
gible role (se e discussion in [20|). In this case, previous 



works 21- 23] have demonstrated an exact decomposition 



of the LDOS as a function of position r and energy E: 
fd 2 K 



p(r,E) 



/7/2-p +°° 



Here, we neglect, for the sake of simplicity, thermal 
smearing effects and the experimental energy resolu- 
tion, which could be implemented straightforwardly. 
A n (R, E) corresponds to the guiding center spectral den- 
sity in LLn, which encodes the complicated dynamics in- 
duced by the disorder. The key object in Eq. (pQ) is the 
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FIG. 2. (color) (a)-(d) Logarithmic absolute values of angular- averaged Fourier transformation (FT) of LDOS obtained in 
the tail of the LLs; thick red lines: experimental data taken from LDOS of Fig. [T^a) (LL0), Fig. Q^c) (LL1), Fig. [^f) (LL2), 
Fig-IHg) (LL3); voltage averaging: 2 mV; thin black lines: FT of numerical simulations for an energy in the tail of the Landau 
band, E — Jtw c • (n — 0.4 + 1/2); the curves are stretched along the y axis in order to match at q = nm -1 ; insets: FT-LDOS; 
upper panels: experiment with dashed half circles indicating q of simulated LDOS minima, lower panels: simulations; (e) 
Angular- averaged FT-LDOS of LL1 with (red curve) and without voltage averaging as derived from Fig.[T]Je) (black), Fig.[TJd) 
(blue), Fig.[T]Jc) (green); (f) Comparison of absolute Fourier-transformed structure factor |F n (q)| (shifted vertically for clarity) 
as given in Eq. (4|) (black lines) and discrete mimima of the curves in (a)-(d); red circles: experiment; blue triangles: simulation; 
black arrow marks |q| = 27r/£t yP = 0.126 nm -1 with £ tyP = 50 nm. 



so-called structure factor 

(-i) n 



F n (R) 



2R/ 



-R 2 /l% 



(2) 



where L n (z) is the Laguerre polynomial of degree n con- 
taining n oscillations (or, equivalently, nodes). 

The structure factor F n (R) is not strictly positive and, 
while associated with the quantum cyclotron motion in 
LLn, cannot be interpreted directly as the wave func- 
tion probability density. Instead, F n (R) corresponds to a 
Wigner distribution [22| , because the physical real-space 
of the guiding center coordinates R = (X, Y) is in fact 
associated to a pair of quantum conjugate variables [24[, 
owing to the commutation relation [X,F] = il 2 B . The 
convolution given by Eq. (pQ) illustrates the difficulty in 
resolving in real-space the nodes of the Landau states 
that are built in the structure factor. Indeed, the sharp 
nodal structure of F n (H) is not only blurred in p(r,E), 
but the smeared nodal patterns must also follow the ran- 
dom meanders of the guiding center trajectories encoded 
in A n (R, E). However, by performing a 2D Fourier trans- 
formation (FT) of the LDOS, a simpler product form is 
achieved 



+oo 

p(q,E) = Y / Fn(<l)A n (q,E), 
where the FT of the structure factor reads: 



(3) 



(4) 



For a smooth disorder potential characterized by a large 
correlation length £ ^> A n (q_, E) possesses disorder- 



induced structures mainly at short wave vectors |q| c± 
27r£ -1 , and should barely vary around the larger momen- 
tum scale |q| ~ R n jl\ oc 1/Z#. Thus, distinguishable 
momentum- variations in the FT of the LDOS (|3]) should 
arise from the universal structure factor (j4]). We now ex- 
amine this issue both using our STS data and numerical 
simulations of a disordered 2DES at high B field. 

In order to reveal the nodal structure of LLs in 
momentum-space encoded into F n (q), we proceed with 
the FT of the LDOS images of Figs, dja) (LL0), Hfc) 
(LLl),[Hf) ( LL2 )> andHIg) (LL3). Since the FT of the 
LDOS displays sign changes due both to the structure 
factor and to the random spectral density in momentum 
space, we focus our discussion on the absolute value of 
the FT signal. We also use a voltage averaging of 2 mV 
of the data in order to improve contrast. The upper 
insets in Figs. EJa)-(d) illustrate the resulting Fourier- 
transformed LDOS (FT-LDOS) images. They exhibit 
a single radial modulation for LL1, changing into dou- 
ble (LL2) and triple (LL3) radial modulations at the 
|q| positions marked by dashed half circles. In order to 
smoothen the random contributions of the spectral den- 
sity, we perform an angular average of the signal as shown 
by the curves (thick red lines) in Figs.[2ja)-(d). A conse- 
quence of this angular averaging is that deviations from 
perfect angular symmetry and noise within the experi- 
ments will add up to a (possibly g-dependent) finite back- 
ground. Therefore the expected zeros from the structure 
factors are shifted upwards and become minima. At LL0 
[Fig. Efa)], the resulting FT curve decreases monotoni- 
cally, composing a disk structure in the FT-LDOS (upper 
inset). However, at LL1 [Fig.[2fb)], the FT curve exhibits 



a dip at |q| = 0.16 nm" 1 followed by a broad hump. This 
results in an additional ring-like structure surrounding a 
smaller disk in the FT-LDOS image (upper inset). The 
number of dips/humps in the FT curve, i.e. the num- 
ber of additional rings in the 2D-image, increases to 2 at 
LL2 [Fig. Etc)] and to 3 at LL3 [Fig.^d)]. This general 
trend is almost quantitatively reproduced by numerical 
Hartree simulations [!, [25| [thin black lines in Fig. Efa)- 
(d) and lower insets] . Note that especially the position of 
the minima in the simulated FT curve shows good agree- 
ment with the position of the dip in the experimental 
one. The simulations diagonalize the wave functions at 
B = 6 T within a random disorder potential calculated 
using TVd and 7V"a of the InSb sample [1, [25| . In order 
to transform the resulting 3D disorder potential into 2D, 
a folding with the confined wave function parallel to B 
deduced from the triangular well approximation [26[ is 
used. The resulting 2D disorder is lower than in the ex- 
periment, since it ignores disorder from the Cs atoms 
[8[. This lower disorder explains the stronger features 
in the simulated FT, however, without any effect on the 
position of the minima. Moreover, the position of the 
minima is robust within each disorder-broadened LL as 
shown in Fig. [2fe). The same dip position is observed 
for all voltages, i.e. for localized as well as for extended 
states. This is remarkable, since the corresponding real- 
space LDOS [Figs.QJc)-(e)] shows very different patterns 
due to the complicated guiding-center motion. The FT 
results confirm experimentally that p from Eq. [3] is ^dom- 
inated by the energy-independent structure factor F n (q) 
at the large momentum scale |q| c± l^ 1 . More precisely, 
the minima appearing at larger q are direct fingerprints 
of the distinct nodal structure within each LL. Fig. [2ff) 
shows a direct comparison of the minima positions de- 
duced from Fig.[2fa)-(d) with |F n (q)|. The mimima from 
the FT of experiment (circles) and simulation (triangles) 
quantitatively match the minima in |F n (q)|. 

We now come back to the real space data, where the 
nodal structure can be observed by appropriate filtering. 
We use the fact that the disorder-induced drift is encoded 
up to |q| ~ 0.13 nm -1 [black arrow in Fig. (2jf)] with our 



potential disorder 27]. The key characteristic feature of 
the real-space nodal structures is encoded at larger |q|. 
Thus, by performing band-pass filtering of the FT-LDOS 
at large q and a subsequent inverse Fourier transform, we 
can identify in real space the transverse nodes decoupled 
from the potential disorder (see Fig. [3]). Since the scales 
of disorder and guiding center motion are not strongly 
different, adequate borders of the filtering are essential to 
improve the visibility of the nodal structure as outlined 
in the supplement [20] . In the numerical simulations, we 
find faint LDOS corrugations around the guiding cen- 
ter trajectories encircling the potential minima, marked 
by crosses in Figs. [3ja) and (b), while sharp oscillations 
perpendicular to the drift trajectory, marked by dotted 
lines, appear for LL1 (two maxima) and for LL2 (three 




FIG. 3. (color online) Real-space nodal structure of LLs along 
the guiding center trajectories; (a), (b) Raw real-space LDOS 
from numerical simulations within tail of LL1 (a) and LL2 (b) 
[E = faj c (n— 0.4+1/2), n = 1,2]. (c), (d) The same real-space 
LDOS as in (a),(b) after Fourier high-pass filtering [28[ within 
|q| > 0.14 nm -1 . Crosses in (a)-(d) mark representative po- 
tential minima, (e), (f) Fourier band-pass filtered real-space 
LDOS from Figs.[T]Je),(f) with additional low-pass filtering |q| 
< 0.52 nm -1 to reduce experimental noise; dashed rectangles 
mark the same area as in Fig. [IJb) ; black dotted lines mark 
double (LL1) and triple (LL2) lines, (g) Line sections along 
the same line marked in (e)(LLl),(f)(LL2) and filtered LDOS 
(LL0) [U from Fig. \V[b) (shifted vertically for clarity). 



maxima) after Fourier filtering [28j [Figs. [3fc),(d)]. The 
double and triple lines are distinct in flat areas of the 
potential (as marked), while in other areas where neigh- 
boring potential maxima are close, interference between 
their corresponding drift states prohibits a clear distinc- 
tion. Within the experimental data, the nodal line struc- 
ture is more difficult to find due to the large number 
of overlapping potential extrema. Careful comparison of 
the data obtained within different LLs in combination 
with Fourier filtering nevertheless reveals some regions 
where the transition from single via double to triple lines 
is observable. For example, consider the marked lines 
(in dashed rectangles) in the filtered LDOS of LL1, LL2 
[Figs.|2e),(f)] and LL0 [Io| [see line sections in Fig.^g)]. 

In summary, FT of spatial STS data leads to a decou- 
pling of the LDOS structures attributed to the quantum 
cyclotron orbit and to the complex disorder-dependent 
guiding-center motion. The former exhibits n minima 
at universal |q| for the nth LL. These minima are asso- 
ciated to the n oscillations of the Landau structure fac- 
tor. By adequately Fourier filtering the real-space LDOS, 
the nodal structure is identified in real space. Our find- 
ings demonstrate that Landau quantization implies dis- 
order independent universal features on the microscopic 
scale. This complements the well known universal quan- 
tum Hall plateaus. 
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